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Abstract 



The effective charge of a rigid polyelectrolyte (PE) approaching an oppo- 
sitely charged surface is studied. The cases of a weak (annealed) and strongly 
charged PE with condensed counterions (such as DNA) are discussed. In 
the most interesting case of the adsorption onto a substrate of low dielectric 
constant (such as a lipid membrane or a mica sheet) the condensed counte- 
rions are not always released as the PE approaches the substrate, because of 
the major importance of the image charge effect. For the adsorption onto 
a surface with freely moving charges, the image charge effect becomes less 
important and full release is often expected. 
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A deep understanding of the adsorption of DNA or other charged biomolecules onto 
oppositely charged membranes is of fundamental importance to understand many key phys- 
iological processes, and many experimental studies have approached this problem from very 
different viewpoints [0], strongly motivated by applications to gene therapy 0. On the 
theoretical side, the adsorption and the interaction of charged rods on a charged surface 
have been studied within the Debye-Huckel approximation in different situations 0. In this 
work, we want to focus on the effective charge of the adsorbed macromolecule. We calculate 
the attraction energy between an infinitely long, charged cylinder (rigid polyelectrolyte or 
PE) parallel to an oppositely charged plane, as a function of their distance h. The en- 
ergy variation leads to the determination of the equilibrium charge density of the rod as a 
function of h. This (effective) charge density can be interpreted in terms of the release of 
condensed counterions for highly charged PE such as DNA, which are beyond the Manning 
condensation threshold, or in terms of the recombination of ionized charges on the rod for 
weak (annealed) PE. We show that in contrast to what would be naively expected, the full 
release of the counterions condensed onto an highly charged rod is not always observed in 
the vicinity of an oppositely charged surface. 

The adsorption energy is derived by perturbating the Gouy-Chapmann solution of the 
Poisson-Boltzmann equation in a planar geometry. The perturbative treatment is strictly 
speaking valid only for low linear charge densities r of the rod Ibt -C 1, where Ib is the 
Bjerrum length Ib = e^/(47re) (all energies are in fc^T units). However, the physical pic- 
ture which emerges from this calculation leads to qualitative statements concerning highly 
charged PE as well. We discuss the adsorption free energy on a substrate of low dielectric 
constant with respect to water {e^ = 80), which is of most practical importance in biology 
related problem (adsorption of DNA on a lipid membrane for which e/p ~ 2 [^]) and other 
situations (adsorption onto the mica surfaces of an SPA: e^c — 6 0). We have checked that 
the case of a membrane of thickness I = 50A and dielectric constant eip = 2 does not show 
quantitative differences with the present situation [l = oo and eip = 0). Finally, we also 
study the case where the charges on the plane are free to adjust to the field created by the 
rod, a situation of great interest for fluid interfaces such as biological lipid membranes. 

The electrostatic potential (j)^^^ near a charged wall of density a > (or equivalently with 
a Gouy-Chapmann length A = 1/{27tIbO')), in a salt solution of average concentration no (or 
Debye length k,~^ with = SttIbUo) satisfies the Poisson-Boltzmann (PB) equation [0. In 
the low salt limit kA -C 1, which is the case discussed in this paper, the potential near the 
wall {kz -C 1) follows the Gouy-Chapmann (G-C) law [0]: 

(1) 

where z is the coordinate normal to the wall (z = at the interface), n^^^ is the density of 
(negative) counterions (or c-i) near the wall. The Gouy-Chapmann solution predicts a dense 
counterion layer (the G-C layer) of thickness A containing a finite fraction of the c-i, followed 
by a diffuse c-i region. The local screening length in the G-C layer is small: ~ \/\/2, 
while the screening length of the diffuse region is self-similar: L,^ ~ z/\/2. For distances 
larger than the Debye length: z > k^^, the electrostatic potential decreases exponentially. 

The potential variation 6(j) due to a negatively charge rod (of linear charge density — r) 
located at an altitude z = h above the plane is calculated by a linear expansion of the G-C 
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theory, provided that the perturbation is smalL A small perturbation of Eq.(|I|): 5n -C n^'^'^ 
(note that 5n is an algebraic average taking into account the sign of the free charges) supposes 
that the perturbation potential 50 is smaller than unity. The calculation of the potential is 
carried out in Fourier space for the coordinate x parallel to the wall and perpendicular to 
the rod: fq = J dxe^'^^f{x). The linearized PB equation and the boundary conditions are: 

and dj(i)\z=h+ - dj(i)\z=h^ = AttIbt (2) 

with natural boundary conditions for z oo 

From the perturbed potential z), the free energy due to the presence of the rod can 
be calculated by a charging process of the rod: 5JF = /J" ((pflf^ + S(f)z=hj dr. The first part 
of the integral gives the interaction energy between the rod and the charged plane and the 
unperturbed G-C counterion layer, while the second part represents the self-energy of the 
rod, and its interaction with the perturbed c-i cloud. The direct interaction energy is: 

5^.„, = 2rlog^^^^ (3) 

The solution of the perturbed PB equation Eq.(||) leads to the "self-energy": 6J-'seif = 
l/2/Br%, with 

rPM (l+p{l + h)y 



(4) 



with p = qX and h = h/X. The integral cutoff is Pm = 2'n-X/a where a is a "microscopic" 
size of order the rod radius (a = 20^4 for DNA P]). This complicated expression can be 
approximated in the two important limits: for /i ^ A 




5-^.e./-i/i.rM| + logU^^ (5) 



and for /i -C A. 




SJ'self ^ -Ibt' -r + + log — (6) 



The self-energy behaves very differently for large and small distances from the wall. It 
shows an attraction between the rod and the plane at large distances (the energy increases 
with the distance), which superimposes to the bare attraction between the two oppositely 
charged macroions. At short distances, there is a strong (logarithmically divergent) repulsion 
between the rod and the plane (Eq.(^)). 

The self-energy shows a deep minimum for a position of the rod which only depends upon 
the G-C length: /imm — 0.8A. This minimum results from a balance between the repulsive 
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image charge effect and tlie attraction due to tlie screening of tlie c-i in tlie Gouy-Cliapman 
layer. Tlie self-energy of a cylinder in a bath of mobile charges is of order /^r^ log L^/ a, where 

is the characteristic ( "screening" ) length of the bath. This energy can be interpreted as 
the interaction free energy of the charges on the cylinder (the interaction having a range 
Lk). Note that although the free ions are mostly of the same sign as the rod, one can speak 
of a "screening" effect, as the interaction between the rod and the unperturbed ion cloud is 
taken into account in 5 Tint (Eq.(^). In addition to the self- interaction, one should consider 
the interaction with the image charge, which in the case e^^o -C e^>o, is a virtual rod of 
same charge located ai z = —h 0. This gives an extra contribution ~ Ibt'^ log Lf^/h per 
unit length. The large and short distance behaviors Eq.(^J^) can now be explained from the 
counterion profile Eq.(|l|), since ~ /i for h > X and ~ A for for h < X. 

The adsorption of a polyelectrolyte onto a biomembrane or other fluid membranes, which 
are generally a mixture of charged and neutral lipid molecules, involves the movement of 
surface charges as a response to the field created by the rod. We address this case, disregard- 
ing the fact that the lipid bilayer is a flexible object which would, to a certain extent, wrap 
around the PE - this interesting phenomenon will be studied in future works. We also 
assume that the charge reorganization at the interface is not limited by the availability of 
moving charges on the plane. The charge on the plane follows a Boltzmann law: a = aQC'^'^ . 
As a result, the boundary condition for the perturbed field on the plane (corresponding to 
Eq.(ID) has to be modified: dz5(t)\z=Q = 47r/Bcro^0|z=o- This affects the self-energy term only, 
which can be expressed similarly to Eq.(^): 6J-'seif = ^I'^lsT'^lp^ with 
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e-'P''- ^ ' ^ - - ' (7) 

\ 3 + 3p + p2 + ) ^ > 

For h > 1 this self-energy is equivalent to the constant surface charge case (Eq.(^), while 
for /i < 1, the approximation gives: 



5-F../->M-r-r^ + iog^ (8) 




The mobile (positive) charges of the surface are attracted toward the rod (x = 0); this 
effectively decreases the "Gouy-Chapman" screening length around the PE and reduces the 
image charge effect. The case of an annealed surface charge is qualitatively similar to (see 
scaling in Eq.(||)), but quantitatively different from the case of a quenched surface charge. 
Since the repulsion of the wall is weakened, the minimum of the self-energy is much deeper, 
and closer to the wall in the case of moving surface charges. 

We now discuss the variation of the line charge of an annealed, or weak, polyelectrolyte 
near a charged wall. The charges on the chain result from a partial ionization of specific 
chemical groups. The ionization occurs at chemical equilibrium with the free ions in solution. 



and is governed by quantities such as the pH of the solution |TT|. Formally, the PE charge 
density can be determined by equating the chemical potential of the charges on the chain, 
to a (given) chemical potential yUo for the free charges. For an infinitely long rigid PE in 
a salt solution, the free energy (per unit length) of the charges on the chain is the sum 
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of the translational entropy of the charges along the rod, the electrostatic energy, and the 
chemical potential: T = r log ra/e + IsT'^^ogna — /igr. The equilibrium charge density 
of the rod for a given chemical potential is obtained by differentiation of the free energy: 
fiQ = log ra + 2lBT\ogKa. The equivalent expression for a charged rod near an oppositely 
charged plane can be computed from the electrostatic free energy Eq.(^,^): /xq = log ra + 
2'rIbt 2p + 21ogfi;(A + h)/2. Because of the minimum in the self-energy, the equilibrium 
charge density r is maximum for a finite height of order the G-C length, and decreases 
sharply near the wall. This non-trivial behavior of the charge density of an annealed PE near 
a wall is mostly due to the importance of the image charge effect in the vicinity of the wall. 
It is extended further below in the case of highly charge PE with Manning condensation. 
Note that the optimum charge of the weak PE can reach higher values if the charges on the 
surfaces are mobile, but still decreases at shorter distances. 

The perturbative treatment is expected to fail in the important case of the release of 
the counterions condensed onto a highly charged rod, as the rod approaches an oppositely 
charged plane. However, the qualitative argument which translates the concentration of free 
charges into a local screening length should still hold in this case. 

A charged cylinder surrounded by its counterions undergoes the so-called Manning con- 
densation [^]. Solutions of the poisson-Boltzmann equation in this geometry |]T3[ predict 
that if the rod is highly charged (namely Ibt > 1), a finite fraction 1 — /3 of the coun- 
terions are confined in the close vicinity of the rod. The electrostatic properties at large 
distances from the rod are the same as those of a rod with an effective charge Ibt* ~ 1. The 
(over) simplified Oosawa picture of counterion condensation |]14| gives a qualitative account 
of this phenomenon [|15|; it is based on a chemical equilibrium between two types of counte- 
rions: condensed c-i with a reduced entropy and subjected to a large electrostatic attraction 
from the rod on the one hand, and free c-i in solution far from the rod on the other hand. 
This picture can be adapted to the case of a cylinder of charge r in a salt solution of density 
uq (or screening length k~^). The effective charge of the cylinder, r* = jSr, is obtained by 
balancing the chemical potentials of the condensed c-i ficond = log 
and of the c-i dispersed among the salt molecules /i/ree = log [^0"^] 
molecule). The resulting fraction of free c-i j3 is: log [81bt{1 — P)] 



{1-I3)tv 



+ 21 BTp log [ko] 

[v is the volume of a c-i 
= il-lBT/3)\og[{Kaf]. 
This result is very similar to the Oosawa relationship for a rod which would occupy a (small) 
volume fraction $ = (/ta)^ in a salt-free solution [jl6|, namely /5 ~ 1 for Ibt < 1 and PIbt — 1 
for Ibt > 1. 

The counterion condensation on a rod near a charged plane can be derived in the same 
way. The chemical potential of the condensed counterions depends upon the electrostatic 
potential on the rod, which can be determined (at the level of the scaling laws) using 
the calculation of the previous sections. Equating the chemical potentials of the free and 
condensed c-i, we obtain the fraction of free c-i (5: 



2lBTp\og 



log [8/Br(l - 
ah 



P)] = log 
+ 2 log 



K{h + X) 



(9) 



where the local screening length follows the asymptotic behaviors: = A for < /i < A and 
= h for X < h < k^^. Examples of the counterion release as a function of the distance 
to the wall is shown on Fig.l (for which proper screening due to the salt has been taken 
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into account for nh > 1). The Manning parameter Ibt = 4 and the radius a = 20A are of 
order those for DNA, three surface charge densities ranging from A = 200A to 50A as been 
considered, corresponding to one charge every (100 A) 2 to {50 Af. The salt concentration 
= lOOOA corresponds to a concentration of 10~^mol/l. While the c-i are released 
as the rod penetrates the Gouy-Chapmann layer, full release {jS = 1) is only reached for 
A = 50y4 - the maximum being 60% for A = 200^4 (it is of the order of 30% for the free rod). 
Furthermore, the free charges recondense at short distance if A is large enough; the effective 
charge reaches 50% of the bare charge for A = 200^4. It can be shown that for large salt 
concentration or weak surface charge of the wall: k^A^ < 2a, the PE in contact with the 
wall can have a lower effective charge than the free PE. 

It should be noted that the short range repulsion due to the image charge can prevent 
real adsorption with an equilibrium contact between the rod and the plane. In many cases, 
however, a short range attraction of non electrostatic origin (such as a hydrophobic force) 
dominates the adsorption. This interaction must be added to the electrostatic free energy 
calculated here in order to study the equilibrium adsorption. 

Different conclusions are reached in the case of moving surface charges. A highly charge 
rod has a strong effect on the surface charge distribution, for the potential it creates is likely 
to dominate over the Gouy-Chapman potential, even near the wall. A quantitative descrip- 
tion of the phenomenon would require the solution of the full non-linear Poisson-Boltzmann 
equation, with complex (non-linear as well) boundary conditions. In the following, we merely 
try to give a feeling of the way moving charges can influence the counterion release of an 
adsorbed polyelectrolyte. We assume that the surface charge distribution obeys Boltzmann 
statistics: ax=o = ctoc""'^'', where the potential (j)h created by the rod at the surface re- 
flects the screening due to the Gouy-Chapmann layer, the image charge, and the fraction 
(1 — P) of condensed counterions: (ph ~ log {L^/hY . Since the local screening length 

Lk — (A + h)/\/2 is influenced by the surface charge a, which in turn, depends upon the 
screening length, we obtain a self-consistent relationship for the local Gouy-Chapman length 
near the rod: 



This expression shows that moving surface charges strongly reduce the screening length near 
the rod, which becomes of order h when the rod is close to the wall (recall that /bt/? ~ 1). 
As a consequence, the interaction with the image charge and the self-interaction along the 
rod are both strongly screened. The dominant interaction between the cylinder and the 
wall is the attractive part given by Eq.(^, and we are likely to observe a full release of the 
condensed counterions in this case. 

To summarize, we have studied the evolution of the effective charge of a polyelectrolyte 
near an oppositely charged plane. The case of a weak PE is studied fairly rigorously, via a 
perturbative treatment of the non-linear Poisson-Boltzmann equation. We show that at large 
distance, the charge of the PE increases as the distance to the wall h decreases, as expected. 
However, the charge decreases "strongly" as the PE enters the Gouy-Chapman layer {h ~ A) 
because of the combined effect of image charge (for the most common case of a wall with a 
low dielectric constant) and self-interaction along the PE. These two effects are very much 
influenced by the value of the Gouy-Chapman length A, and are partly suppressed in the 




(10) 



6 



case of a fluid interface with moving surface charges (for a fluid lipid bilayer for instance), 
where the charge of the adsorbed PE can reach higher values. 

The most interesting case of a strongly charge PE beyond the Manning condensation 
threshold (such as DNA) is discussed qualitatively, using "scahng" arguments inferred from 
the perturbation theory. We predict that in the case of fixed surface charges, and in contrary 
to a widely spread idea, most of the condensed counterions are not released if the Gouy- 
Chapman length is larger than the radius of the rod: A ^ a. In the case of freely moving 
surface charges, a full release of the condensed counterions is expected, as the effective 
Gouy-Chapman length near the rod is of order the rod radius. 

In all the discussion, we have assumed that the Gouy Chapman length is larger than a 
molecular size. In many real cases, the two lengths are of the same order of magnitude and 
the finite size of the ions must be taken into account in order to obtain quantitative results. 
In this case our results can at best be considered as qualitative. 

We would like to thank A. Johner and R. Netz for many stimulating discussions. 
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FIGURES 




FIG. 1. Fig.l Fraction (3 of free counterions as a function of the distance to the wall in A 
units (dashed: (3{h — > co)) for a = 20A, Ibt = 4, = lOOOA, for decreasing values of the 
Gouy-Chapmann length: A = 200^4, 100^ and 50^ 
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